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Abstract. We revisit the single item continuous-time inventory model. In the same problem setting as 
, in Bensoussan et al. [6], we show the optimality of an (s, <S')-policy for a general spectrally positive Levy 

. demand process. Using the fluctuation theory of spectrally one-sided Levy processes, we express the value 

^Sj . function analytically using the scale function. Numerical examples under a Levy process in the /3-family 

^ I with jumps of infinite activity are provided to confirm the analytical results. The case with no fixed ordering 

cost is also studied. 
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1. Introduction 

In this paper, we study the inventory control problem where the objective is to minimize the sum 
^ . of the surplus cost incurred over time and the ordering cost incurred at each replenishment. The latter 
&^ '. consists of the fixed and proportional costs. Existing literatures suggest that an optimal ordering policy 
^ is in many cases of (s, 5')-type; namely, it is optimal to bring the surplus level up to S whenever it goes 
■ below s. This is a well-studied problem with many variants. Typically, the underlying stochastic process 
2 ■ is modeled using a compound Poisson process and/or a continuous process. In particular, in a seminal 
^ ■ paper by [6], they studied the case where the demand is a mixture of a Wiener process and a general 
compound Poisson process. The optimality of the (s, S')-policy in their problem setting was proved with 
a complement by [4]. See [6] and [25] and references therein for a detailed review of inventory control 



^ _ problems. 



This paper studies a generalization of the existing literature where the demand is a general spectrally 
positive Levy process, or a Levy process with only positive jumps that is not a subordinator. This 
generalization encompasses jumps of infinite activity and/or infinite variation that had not been studied 
previously in inventory models. It accommodates important processes such as the variance gamma, 
CGMY and normal inverse Gaussian processes as well as classical examples such as the gamma process 
and a subset of stable processes. The ones perturbed by Wiener processes are also accommodated. 
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As opposed to the quasi- variational inequalities (QVI) approach as in [6], we take advantage of the 
fluctuation theory of spectrally one-sided Levy processes (see e.g. [7, 9, 17]). Using the so-called scale 
function, the expected total costs corresponding to the (s, S')-policy for s < S can be written explicitly. 
The values of s and S are then determined by the continuous (resp. smooth) fit at s when the demand 
process is of bounded (resp. unbounded) variation together with another usual condition of the slope 
at S. The optimality is verified using the continuous/smooth fit at s, the martingale and some known 
properties of the scale function, and also the results by [4]. The optimal levels (s*, S*) are the zero of the 
function Q(s, S) define in (3.16) below and that of its derivative with respect to the second argument. The 
resulting value function can be written concisely via the scale function, which has analytical expressions 
for many cases (see, among others, [14, 17]) and in general can be approximated efficiently by, e.g.. 



We also consider the case without a fixed ordering cost, or the case with solely a proportional ordering 
cost. We shall show, in this case, the optimal policy is of barrier type. Using the fluctuation theory 
of reflected Levy processes [1, 24], the value function can again be written concisely using the scale 
function. 

In order to confirm the analytical results obtained in this paper, we give numerical examples with a 
quadratic surplus cost and a demand process in the /3-family of [15] with jumps of infinite activity. We 
see that the optimal levels (s*, S*) and the value function can be computed instantaneously. We analyze 
the results and also confirm the convergence as the fixed ordering cost decreases to zero. 

To our best knowledge, the inventory control problem has not been solved for a general spectrally 
positive Levy demand process in this generality. As related literature, the optimal dividend problem 
under transaction costs is solved for spectrally negative and positive Levy processes by [20] and [3], 
respectively. Their techniques are similar to ours; the optimality of the (s, S')-policy in their settings 
are shown via the fluctuation theory and the scale function. For related singular control problems, see 
[1, 19]. For a detailed review on the inventory models and impulse control, see [5,21, 22]. 

The rest of the paper is organized as follows. Section 2 gives a mathematical model of the problem. In 
Section 3, we compute the expected total costs under the (s, 5')-policy via the scale function. Section 4 
shows via the continuous/smooth fit principle to obtain a candidate policy and show its existence. Section 
5 verifies its optimality. Section 6 studies the case without a fixed ordering cost. We conclude the paper 
with numerical results in Section 7. 

2. Mathematical Formulation 

Let (fi, J^, P) be a probability space hosting a spectrally negative Levy process X = {Xt;t > 0} 
whose Laplace exponent is given by 



[11,26]. 



(2.1) 




(0,oo) 



—sz 
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where z/ is a Levy measure with the support (0, oo) that satisfies the integrability condition J^^ ^^(1 A 



(dz) < oo. It has paths of bounded variation if and only if cr = and f^^ z iy{dz) < oo; in this 
case, we write (2.1) as 

iP{s) = Ss+ [ {e-'' -l)iy{dz), seR 

J{0,od) 

with 5 := 7 + zu^dz). We exclude the case in which X is a subordinator (i.e., X has monotone 
paths a.s.). This assumption implies that 5 > when X is of bounded variation. Let be the conditional 
probability under which Xq = x (also let P = Fq), and let F := {J't : t > 0} be the filtration generated 
byX. 

We model the demand process as its dual —X, or a general spectrally positive Levy process. We 
assume that Xi has a finite moment. 

Assumption 2.1. We assume that /i G (— oo, oo) where 

(2.2) fi:=E[X,]=ij\0+). 

An (ordering) policy n := {LJ; t > 0} is given in the form of an impulse control (T^, u'^; T2, ^2 ; ■ ■ ■ ) 
with = Yl,rT''<t'^l where {T^; i > 1} is an increasing sequence of F-stopping times and > is 

■ i — 

an J^Ti-i^s^s arable random variable for i > 1. Corresponding to every policy tt, the surplus process is 
given by U'' = {U]^;t >0} where 

■= Xt + L^, t> 0. 
The problem is to compute, for g > 0, the total expected cost given by 

/•oo ^ 

/ e-^V(f/r)dt + 5^e-^^^^7«) 
^>^o i=l ^ 

POO 

/ e-'^'fiUndt+ J2 e-''*^?(ALni{ALr>c 



w^(x) := Er, 



X G 



with AL^ = LJ — LJl, and to obtain an admissible policy that minimizes it, if such a policy exists. 

Here, / : M t-> [0, 00) corresponds to the cost of holding and shortage when x > and a; < 0, 
respectively. Regarding g, we assume 

(2.3) g{y):=Cy + K, y > 0, 

for some unit cost of the item C > and fixed ordering cost K > 0. We shall study the case K = 
separately in Section 6. As in [4, 6], we assume the following. 

Assumption 2.2. (1) / is a piecewise continuously differentiable function with /(O) = 0. 
(2) There exists a number a such that the function 

(2.4) f{x) := fix) + Cqx, x G M, 
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is increasing on (a, oo) and decreasing and convex on (—00, a). 
(3) There exist a cq > and an xo> a such that f'{x) > cofor x > xq. 

As in [4, 6], this is a crucial assumption for our analysis. This implies Lemma 4.2 below that will be 
used to verify the existence of the optimal policy. 
Finally, the (optimal) value function is written as 

(2.5) v{x) := inf Vn{x), x G M, 

ttGII 

where n is the set of all admissible policies. 

3. The (s,S)-policy 

We aim to prove that the (s*, S'*)-policy is optimal for some —00 < s* < 5* < 00. For —00 < s < 

5 < 00, an (s, S')-policy, Hs^s ■= {Lt'^', ^ > O}, brings the level of the surplus process If^'^ := X + L^'^ 
up to S whenever it goes below s. Let us define the corresponding expected total costs as 



(3.1) v,,{x) := E, 



POO 



'>0} 

0<t<oo 



In this section, we shall express these in terms of the scale function. 

3.1. Scale functions. Fix q > 0. For any spectrally negative Levy process, there exists a function called 
the q-scale function 

W^"^ : M ^ [0,00), 

which is zero on (— oo,0), continuous and strictly increasing on [0, 00), and is characterized by the 
Laplace transform: 

/•oo 1 

/ e-'''W^'^\x)dx = — , s > <^(q), 

Jo ^(s) - q 

where 

$(g) := sup{A > : tpiX) = q}. 

Here, the Laplace exponent ^/^ in (2.1) is known to be zero at the origin and convex on [0, 00); therefore 
$(g) is well-defined and is strictly positive as g > 0. We also define, for x G M, 

W^''\x) := [ W^''\y)dy, 
Jo 

Z^i^x) := l + qW^''\x), 

Z^'^x) := / Z^'>\z)dz = X + q / W^''\w)dwdz. 
Jo Jo Jo 

Because W^"'^ is uniformly zero on the negative half line, we have 

(3.2) Z^''\x) = l and Z^''\x) = x, x < 0. 
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Let US define the first down- and up-crossing times, respectively, of X by 
(3.3) := inf {t > : Xt < b} and r+ := inf {t > : > 6} , be R. 

Then, for any 6 > and x < b. 



(3.4) 



1 



E, 



Z^''\b) 



X 



W(i\b) 



X) 



Fix A > and define V'a(-) as the Laplace exponent of X under P'^ with the change of measure 

exp{\Xt-'4j{\)t), t>0; 



dP 



see page 213 of [17]. Suppose W^''^ and Z]^'''' are the scale functions associated with X under P^ (or 



equivalently with ^x{-)). Then, by Lemma 8.4 of [17], 



X) 



e~^^iy('')(x), a; G M, which is 



well-defined even for q < ^/'(A) by Lemmas 8.3 and 8.5 of [17]. In particular, we define 



W^^g){x) 

which is known to be an increasing function. 



X) 



X G 



Remark 3.1. (1) If X is of unbounded variation or the Levy measure is atomless, it is known that 
W(^) is C\R\{0}); see, e.g., [8]. Hence, 

(a) Z^^) is C^(R\{0}) and C°(]R)/or the bounded variation case, while it is C^(M\{0}) and 

(R) for the unbounded variation case, and 

(b) Z*'^^ is C^(M\{0}) and C^{R) for the bounded variation case, while it is C^(M\{0}) and 

(R) for the unbounded variation case. 
(2) Regarding the asymptotic behavior near zero, we have that 

, if X is of unbounded variation, 
|, if X is of bounded variation, 

if a > 

(3.6) 



(3.5) 



^^^(o) 



M/(9)'(0+) := limiy^")'! 



X) 



oo, 

q+u{0,OD) 



(3) As in (8.18) and Lemma 8.2 o/[17], 



ty(<?)(y) - 1^(9) (x) 

In all cases, iy('?)'(x-) > W^'^y {x+) for all x G 



if a = and z/(0, oo) = oo, 
ifX is compound Poisson. 



y > X > 0. 
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As necessary tools for subsequent discussions, we shall define, for all x > 0, 



(3.7) 



e('?)( 



x] 



X- 



9)ra;) = e*('')^iy4(,)(x+)>0, 



e^'^a:) := W'^''\x) - ^q)W^''\x) > 0. 

Here the positivity of the former holds because — ^ x > and that of the latter holds 

because it is an integral of the former. Probabilistically, the former is known as the g-resolvent measure of 
the descending ladder height process of X (see, e.g., [18]); their positivity will be important in deriving 
the existence of the optimal solution and the verification of optimality. 

3.2. The expected total costs under the (s, S')-poUcy. Now we obtain (3.1) using the scale function. 
By (2.3) and the strong Markov property, it must satisfy, for every S > s and x > s. 



(3.8) Vs,s{x) = E, 

Following [4, 6], we define 

(3.9) 

With 



e-'"f{Xt)dt 



+ E, 



-qrs 



Vs,siS). 



Vs,six) := Vs,six) + Cx, X e 



-qt 



f{Xt)dt 



(3.10) k{s,x):=E^ 
and using (3.4), we can write (3.8) as 

(3.11) Vs,six) = k{s,x)+E^ 
where 

(3.12) Vs,s ■= Vs,s{S) 



1 



+ KE^ 



Cx^ X ^ 5, 



Vs,s = k{s,x) - 

1 Q^''\S-s) 



$(g) 



X > s, 



S> s. 



Here the last equality holds by solving (3.11) for x = S. Hence our first task is to compute /c(-, ■); 
together with (3.12) we can then obtain (3.11). 

The first expectation of (3.10) can be obtained using the g-resolvent measure written in terms of the 
scale function; see, e.g., [12, 17]. Define, for any measurable function h and s E M, 



(3.13) 



vl>(s;/i) := 

iPs{x; h) :-- 



-HQ)y 



h{y + s)([y 



'HQ){y-s) 



I W^'i\x-y)h{y)dy, xe 

J s 



Here ips{x; h) = for any x < s because W^'^^ is uniformly zero on (— oo, 0). Note that \E'(-; /') is 
equivalent to (4.23) of [6] (times a positive constant). 

We assume the following, which is an assumption weaker than the polynomial growth assumption of 
[6]. This guarantees that the expectation in Lemma 3.1 below is well-defined and finite. 
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Assumption 3.1. We assume \l/(s; /) < oo for s G M. 

Lemma 3.1. We have [ e-«*/(Xt)dt] = W^'i\x - s)*(s; /) - j) for any x, s G M. 

Using the following lemma, we shall write it in terms of / as in (2.4), which will make our analysis 
easier. 



Lemma 3.2. For s G 



Cq 



Vs{x- f) = <^,(x; /) - C sZ^''\x -s) + Z'^'\x - s) 



X 



X G 



Proof. The first claim holds by integration by parts. For the second claim, we have 

^s{x; f) - ^s{x- /) = -Cq [ yW'^''\x - y)dy 



-Cq 



sW^'\x-s) 



-C 



-X 



+ sZ^''\x - s) + Z^''\x - s 



□ 



Using Lemmas 3.1 and 3.2, for x, s G M, 
(3.14) 



e~'^'f{Xt)dt\ =e^''\x-s) 



Cq 



+ ^q)W'^\x-s 



$(g) V$(g) 
Cq 



^{qy 



Z^^\x — s) — (x — s) 



7il), 



Regarding the second expectation of (3.10), Proposition 2 of [1] gives Ea;[e '^^o X^^] = Z (x) 
fjW^''\x) - (g - n^{q))W^''\x)/{^{qf). This together with (3.4) gives the following. 



Lemma 3.3. For x, s G 



E, 



7il), 



[X — s] — s 



/i\ q 



Q^'\x-s) 



q 



.W(i\x-s). 



q J <l>(g) $(g)2 
Substituting (3.4) , (3.14) and the term in Lemma 3.3 in (3.10), we can write k(-,-) as follows. 

Proposition 3.1. For any x > s, 



(3.15) k{s,x) = e^''\x-s) 

where we define 



q 



$(g) 



q 



Q{s,x) 



(3.16) 



7il) 



g{s, x) := $(g)*(s; f)W'{x - s) + K - y.,(x; /), x > s. 
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4. Candidate policies 

Using the results in the previous section, we now have an analytical expression for (3.11). By (3.12) 
and (3.15), 

(4.1) J^^s,^J^^ns;f)-J^(K^^)^ 



Substituting this and (3.15) in (3.1 1), 

q(9)/^ — s) 

(4.2) Vs,s{x) = ~^7-^ — —g{s,S) + g{s,x)+Vs,s, x>s. 

4.1. Continuous/smooth fit. In order to narrow down the candidates of {s,S), we use the continu- 
ous/smooth fit principle at s. This together with an additional condition at S will give a pair of the 
candidate policy. It is well-known, in the existence of a diffusion component, that smooth fit holds for 
impulse control problems; see [13]. On the other hand, for a Levy process of bounded variation, contin- 
uous fit may be used alternatively. This is well-studied particularly for optimal stopping problems; see, 
e.g., [10, 18]. Here we apply continuous fit for the case X is of bounded variation and smooth fit for the 
case it is of unbounded variation. 
By taking x I s in (4.2) and because 

(4.3) g{s,s+) = K >0, sgM, 
we have 

vs,sis+) = ^ ^ Gis, S) + K + vs,s = ^ ' g{s, S) + €s,s{s-)- 

Note that 6^''''(0+) = if and only if X is of unbounded variation in view of Remark 3.1(2). Hence, 
the continuity at x = s holds if and only if ^(s, S") = for the case of bounded variation while it holds 
automatically for the unbounded variation case. 

For the case of unbounded variation, we further pursue the differentiability at x = s. Differentiating 
(4.2), 

(4.4) (x) = -—-^ ^-g{s,S) + n{s,x), x>s, 

where, for all x > s, 

d 

'H{s,x) := —g{s,x) 

(4.5) ox 

= ^q)^{s; f)W^''\x -s)- ^'^{x; /) = ^{s; f')W^'^\x - s) - <^,(a;; /'). 
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Here, the last equality holds because by integration by parts, 




and ^{y; f) = [f{y) + ^{y; /')]/$(?)• Because 'H{s, s+) = for the case of unbounded variation, we 



see that the differentiability holds if and only if ^(s, S') = as well. 

We now turn our attention to the slope at S. As a natural condition that is used typically in impulse 
control problems, we pursue a candidate for which the slope of Vs,si-) ^ zero. If we impose 
g{s, 5) = in (4.4), we have v'^^siS) = ^(s, S). Hence ii^siS) = if and only if S) = 0. 

We summarize the results in the lemma below. 

Lemma 4.1. Suppose (s, S) are such thatQ{s, S) = Tii^s, S) = 0. Then 

(1) Vs^s is continuous (resp. differentiable) at s when X is of bounded (resp. unbounded) variation, 

(2) vi^siS) = 0. 

4.2. Existenceof (s*,^*). We shall show that there exists a pair (s*, S"*) such that ^(s*, S**) = n{s*,S*) = 
0; in the next section we show that this requirement is sufficient to prove the optimality of the (s*, S*)- 
policy. 

Toward this end, we shall use the following results by [6], which hold directly from Assumption 2.2. 

Lemma 4.2 (Proposition 5.1 of [6]). (1) There exists a unique number qq < a such that ^'(ao; /') = 
0, '^{x; /') <Oifx< flo and /') > Oifx > oq. 
(2) ^'(x;/') >Oforx<a. 



(3) ^{x-f ) > co/<l>{q)forx>xo. 



In order to use this lemma efficiently, we first rewrite (3.16) and (4.5) as integrals of /'). 



Lemma 4.3. For x > s, 



'H{s,x) = ^{x-f')W^''\0)+ / ^{yj')e^'\x-y)dy. 




Proof. We shall show the former because the latter is simply a derivative of the former. 
By integration by parts and because "^{y; f) = [f{y) + '^{y; f')]/^{q), 

vl>(s; f)W^'\x -s)= r hiy; f)W^^\x - y) - vl>(y; f'W^'\x - y)] dy 





$(g) 
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Substituting this in (3.16), we have the claim. 

Lemma 4.4. (1) For every fixed s G M, lim^-t-oo G{.s, S) = oo. 
(2) For every fixed S* G M, lim^^oo G{s, S) = —oo. 

Proof. (1) By Lemmas 4.2(3) and 4.3, for any S* > s V Xq, 

G{s. ^) > ^ £ - y)^y + f'^" /')e^'^(^ - y)dy + k. 

The claim is now immediate because, for e > and sufficiently large S, 

e^^V - y)dy = / e^'\y)dy > / e^'\e)dy ^ oo, 

_Va,'o Jo Jt 

where the inequality holds because 9*^'^'' is positive and monotonically increasing. 
(2) Similarly, if we choose h < oq, we have, for any s < S Ab, 

gis,S)= r'^iyJ')e^'\S-y)dy+ f ^{y- ~r)Q^'\s - y)dy + K. 

Js JSAb 

By Lemma 4.2(1) and (2) and because /') < 0, we have, for any e > and sufficiently small s, 

fSAb , , rShb 



□ 



/ ^{y; nQ^'\s - y)dy < ^{b; f) / Q^'\s - y)dy 

J s J s 



= ^{b; h f ' Q^'\y)dy < ^{b- f) f ' Q^'\S - 5 A 6 + e)dy ^ -oo. 

JS-SAb Js-SAb+e 

□ 

We are now ready to show the existence of (s*, S*). While this is shown analytically shortly below, 
numerical plots of Q and "H in Figure 1 of Section 7 are certainly helpful in understanding the proof. 

For any S > s > uq. Lemmas 4.2(1) and 4.3 imply that Tiis, S") > uniformly. This together with 
(4.3) implies that Gi-s,-) starts at K > (at 5 = s) and increases monotonically as S* t C)0 while never 
touching the zero line. 

Let us now start at s = and consider decreasing the value of s toward — oo. By Lemma 4.4(1) 
and because 'H{s,S) < for any S < a^, S{s) G mis>sGis, S) > ao always exists and attains a 
local minimum (hence 7/(5, S'(s)) = 0). Moreover, Lemma 4.4(2) implies, for sufficiently small s, 
inf5>s G{s, S) < 0. Hence, there must exist 



(4.6) s* := sup <^ s < ao : inf g{s,S) = 0} and S* G arg inf G{s, S), 

I S>s J S>s 

which clearly satisfy G{s*,S*) = 'H{s*,S*) = 0. Note that while s* is unique by construction, S* 
may not be unique. Moreover, by construction, s* < ao and S* > (because ^(s*, S) is negative for 
S G {s*, ao)). This completes the proof of the following. 
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Proposition 4.1. There exists s* < a^andS* > ao such that (4.6) holds with His* , S*) = G{s*,S*) = 0. 

5. Verification of optimality 

With the (s*, S*) obtained in Proposition 4.1, we shall show that the (s*, S'*)-policy is optimal. By 
substituting g{s*, S*) = in (4.2), 

(5.1) Vs*,S*{x) =G{S*,X) +Vs*,S', x>s*. 

Notice from (4.1) that 

q q 
Substituting this and by the definition of ■) as in (3.16), 

(5.2) is*,S'{x) = ^^{s*- j)Z^'i\x - s*) - i^Ax- f) - ^, x> s*. 

q q 

The function can be recovered by (3.9). By Lemma 3.2, we can also write 

(5.3) Vs^M^) = (^*(^*; /) + ^) Z^'\^ -sn-C (z^'\x - s*) + - /). 

It is clear from (3.2) that this expression is also valid for x < s*. 

Let C be the infinitesimal generator associated with the process X applied to a sufficiently smooth 
function h 

Ch{x) := -fh'{x) + -a'^h"{x) + / [h{x - z) ~ h{x) + h'{x)zl{o<z<i}] f^idz), x eM. 
2 Jo 

By Remark 3.1(1) and Lemma 4.1, the function Vs*,s* is C°(M) and Ci(M\{s*}) (resp. C\R) and 
C^(]R\{s*})) when X is of bounded (resp. unbounded) variation. Hence, Cvs*,s*i-) makes sense any- 
where on R\{s*}. 

Proposition 5.1. (1) (£ — q)vs*,s*{x) + f{x) = Ofor x > s*, 
(2) (£ - q)vs*,s4x) + fix) > Oforx< s*. 

Proof. (1) By Proposition 2 of [1] and as in the proof of Theorem 8.10 of [17], the processes 

^-.(*A.,A.+)^(.)(^^^ ^ and e-''(*^vA^^)i?(.)(x, ,*), ^ > q, 

S -t> s -t> 

for any B > s* and R^'^\x) := Z^''\x) + fi/q, x eR, are martingales. Thanks to the smoothness of Z*^"^^ 
and Z^^^ on (0, oo) as in Remark 3.1(1), we obtain (£ - q)R^'^\y - s*) = [C - q)Z^'^\y - s*) = for 
any y > s*. This step is similar to the proof of Theorem 2.1 in [2]. 
On the other hand, as in the proof of Lemma 4.5 of [12], 

{C- q)^s'{x]f) = f{x). 

Hence in view of (5.3), (1) is proved. 
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(2) For X < s*, because Cvs*,s* is a constant and by (3.2) and (5.3), we have 

(5.4) ^ [(£ - g)tv,5. (x) + fix)] =qC + f\x) = f\x) < 0, a.e. x<s*, 

where the last inequality holds because x < s* < ciq < a and by the definition of a. Hence {C — 
1)^3*, s*{x) + f{x) is decreasing on (—00, s*). Therefore it is sufficient to show that 

(5.5) iC-q)vs*,s*is*-) + f{s*)>0. 

Suppose X is of bounded variation. Because s* < (and hence /') < 0), we have 

{5^, ,..(s*+) = nis*, s*+) = /')1^(«)(0) < 0. 

This implies, together with (1), continuous fit as in Lemma 4.1 and 6 > 0, shows (5.5) for the case of 
bounded variation. 

Suppose X is of unbounded variation. Then 

<,,5.(x) = n'is^x) = ^ £ ^{y; /')e('?)(x - y)dy, x > s*. 

Taking limits and notice again that s* < ao, 

<.,5.(s*+) = vi>(,-/')e(«)(o+)<o. 

This implies together with (1) and smooth fit as in Lemma 4.1 shows (5.5). □ 
Proposition 5.2. For every x eM., we have fs*,s'*(x) < K + inf„>o [Cu + Vs*,s*{x + '^^)]- 

It is clear that showing the above is equivalent to showing Vs*^s* (x) < K + inf„>o ^s*,^* {x + u). 

Lemma 5.1. The following holds true. 

(1) Vs',s'-iS*) = inf^eK{;,.,s.(x). 

(2) Vs*,s*{x) is decreasing on [s*,ao]. 

(3) lim^^ooV s*,s*{x) = 00. 

Proof. By (5.1) and because S* minimizes G{s*, x) over x as in Proposition 4.1, the first claim holds. 
The second claim holds because f/. 5, (x) = 'H{s*,x) < on [s*, ao] in view of Lemmas 4.2(1) and 4.3. 
The third claim holds by Lemma 4.4. □ 

Forx < s*. Proposition 5.2 holds immediately by Lemma 5.1(1) and because f;^,. 5. (x) = Vs*,s'-{s*) = 
Vs*,s*iS*) + K. For s* < X < ao, it also holds by Lemma 5.1(2). The proof for x > ao is the most 
difficult part of the verification. However, thanks to Proposition 5.1 and Lemma 5.1, we can apply 
Theorem 1 of [4] where they used the concept of the non-i^-decreasing function. While they assume a 
mixture of a Wiener process and a compound Poisson process, the arguments are still valid for a general 
spectrally negative Levy process. Hence, we have Proposition 5.2. 
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By Propositions 5.1 and 5.2, we have shown that Vs*,s* satisfies the QVI. The rest of the proof of 
optimality is omitted because this is a well-studied problem and the results hold by applying Ito's formula 
to Vs*,s*iU^) for an arbitrary policy vr G H. Here one needs to apply an appropriate version of Ito's 
formula (e.g. Theorem 3.2 of [23]) since the value function Vs*^s* is not smooth enough at s* to apply 
the usual version. However, thanks to the continuous/smooth fit as in Lemma 4.1, the local time terms 
that might accumulate around s* will vanish and hence have no effects after all. See also Theorem 6.2 of 
[22] for the verification lemma for impulse control problems. 



Theorem 5.1. The (s*, S*) -policy is optimal and the value function is given by Vg 



as in (5.3). 



6. The case with no fixed ordering cost 

In this section, we study the variant of the problem with K = 0. Here, we widen the set of policies to 
accommodate also the processes containing diffuse components; we consider the set of vr := {L^; t >0} 
given by a nondecreasing, right-continuous and F-adapted process starting at zero. With U[ := Xt + LJ, 
t >0, the problem is to compute the total costs: 



E.,. 



e-'^'ifiUndt + CdL^ 



X G 



for some C > and to obtain an admissible policy that minimizes it, if such a policy exists. 

From the results in the previous section, the optimal policy is easily conjectured. We have seen, for 
the case K > 0, an (s*, S'*)-policy is optimal for some s* < < S*. Because the distance between s* 
and S* is expected to shrink as K decreases, it is a reasonable guess for the case K = that a barrier 
policy with lower barrier oq is optimal. We shall show that it is indeed so. 

Define, for s G M, 



LI 



sup [s 

0<t'<t 



Xh) V 0, t> 0. 



The corresponding surplus process := Xt + becomes a reflected Levy process that always stays 



at or above s. Our objective is to show that Va 



infjren ^^.{x) for all x G M where 



(6.1) 



VAX 



e-"* ifiUndt + CdLt) 



x, s G 



and n is the set of all admissible policies. 

The fluctuation theory of the reflected Levy process has been well-studied. The expression (6.1) can 
be computed easily via the scale function. 



Lemma 6.1. We have E,^ 



Jo e 



-Z^'\x 



s) — jJ'/q + Z'^''\x — s)/^{q)for any x,s E 



Proof. As in the proof of Theorem 1 of [1], if we define := inf {t > : f/° > 5} for 5 G 



E, 



Z(i){B) 
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with/(x) 



(X 



FromExercise8.5of [17]andl'H6pital'srule, Z^'^\b)/Z^'}\B) 



— > 



$(g) ^ as 5 t oo. Hence, by the monotone convergence theorem, 



JO 



-l{x) + Z^'i\x) lim 



-l{x). 



By shifting the initial position of X, the proof is complete. □ 
Lemma 6.2. We have E^. [ e-^^f{Uf)dt] = Z'^'^^x - /)/g - </?,(x; f) for any x,seR. 

Proof. By Theorem l(i) of [24], for every B > x 



E,. 



e-'?7([//)dt 



(s,oo) 



(s,oo) 



Z^'\x - s)^^^^^^ - W^'^\x - y)] f{y)dy 

Z^^Hx - .)SE^e-(^)(-)?^M^^ - Wi''\x - y)]fiy)dy. 



Thanks to Assumption 3.1, the boundedness of W^'^\-) / Z^'i\-) in view of (3.4) and because W^(^q) is 
increasing, this is bounded uniformly in B. Hence, the dominated convergence theorem along with the 
convergence W^'^^B) /Z^'^\B) -)■ /q (by Exercise 8.5 of [17]) yields the resuk. □ 



Combining the two lemmas above, we can now write 
(6.2) vaoix) = -C (^^'\x - ao) + ^) + ' «o) 



^(ao;/) ' 



$(g) 



which holds also for X < sby(3.2). By Lemma 3.2 and because \l/(ao; /) = /(ao)/$(g) by the definition 
of ao, we can write 

(6.3) Vaoix) := Vaoix) + Cx = + — — — —f{ao)-(pao{x;f), xeR. 

It turns out that our choice of s = oq guarantees the smoothness of (6.2) (or equivalently (6.3)) at 
that is even stronger than the case K > as in Section 4.1. First, it is clearly continuous. Taking a 
derivative, because (p'a^ix; f) = W^''\x - ao)/(ao) + ^ao{x; /'), 

(6.4) <„(x) = W^^\x - ao)/(ao) - /) = /') ^ 0. 
In particular, for the unbounded variation case. 



KXx) 



W^''y{x-y)f'{y)dy 



> 0. 



ao 



These results together with the smoothness of the scale function as in Remark 3.1(1) show the following. 

Lemma 6.3. The function is C^(]R) (resp. C^{'R))for the case X is of bounded (resp. unbounded) 
variation. 
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We shall now prove that solves the variational inequality. By Lemma 6.3, Cva^ix) makes sense 
for all X G M. 

Proposition 6.1. (1) {C — q)vao{x) + f{x) = Ofor x > ao, 
(2) (£ - q)va,{x) + fix) > Oforx < ao- 

Proof. (1) In view of (6.2), the proof is similar to Proposition 5.1. 

(2) Using X < ao < a and the definition of a as in Assumption 2.2, the same steps as in the proof 
of Proposition 5.1(2) yield that (C — q)vaf){x) + f{x) is decreasing on (— oo, ao). Now the claim holds 
immediately by(l)andLemma6. 3 because (£—g)t'ao(ao—) + /(ao) = {^—Q)vao{o,o+)+ f (0,0) = 0. □ 

Proposition 6.2. We have v'^^^ (x) = Ofor every x < ao and v'^^ (x) > Ofor every x > ao. 

Proof. Recall (6.4). When x < ao, it is clear that (pao{x; /') = 0. When ao < x < a, because f'{y) < 
for every ao < y < a, the claim holds in view of the definition of /^aoi^] /') as in (3.13). 
Now suppose X > a. We write 

(6.5) /') = e*^'')^ f e-*('')W$(,)(x - y)f'{y)dy. 

J ao 

We shall show that the integral is negative. Because f'{y) > for any y > x > a, 

n ^ poo 

(6.6) / e-*(«)V~'(2/)d2/ < / 6""^'^^^ f'{y)dy = e-''^'^^^°^{aoJ') = 0. 

J ao J ao 

Because W<s,(q) {x — y) is decreasing in y, the inequality of (6.6) implies that the integral in (6.5) must be 
negative. Indeed, in the integral, W^(^q){x — y) puts less weight for y E [a, x) where f'{y) > and more 
weight for y G [oq, a) where f'{y) < 0. This completes the proof. □ 

By Propositions 6.1 and 6.2 and thanks to the smoothness as in Lemma 6.3, a standard verification 
lemma (see e.g. Proposition 4 of [1]) shows the optimality. 

Theorem 6.1. The barrier strategy is optimal and the value function is given by (■) as in (6.2). 

7. Numerical Results 

In this section, we conduct numerical experiments using the spectrally negative Levy process in the 
/3-family introduced by [15]. The following definition is due to Definition 4 of [15]. 

Definition 7.1. A spectrally negative Levy process is said to be in the P-family if (2.1) is written 

i^{z) = 5z+ \a^z' + ^ \B{a + |, 1 - A) - B{a, 1 - A)| 

for some 5 G M, a > 0, /3 > 0, > 0, A G (0,3)\{1,2} and the beta function B{x,y) := 
T{x)T{y)/T{x + y). 
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This is a subclass of the meromorphic Levy process [16] and hence the Levy measure can be written 

oo 

i=i 

for some {pfc, 77^; k > 1}. The equation 'ip{-) = q has a countable negative real-valued roots {—^k,q', k > 
1} that satisfy the interlacing condition: ■ ■ ■ < —rjk < —^k,q < ■ ■ ■ < — 7/2 < —■^2,9 < —Vi < ~Ci,<j < 0. 
By using a similar argument as in [1 1], the scale function can be written as 

(7.1) W^^^Ha:) = ^^;^-E^-^^^^^^^^^ -^0' 

where 

q *(?) + fi,, V 1i / iil-r^ 

We suppose 6 = 0.1, A = 1.5, a = 3, /3 = 1 and w = 0.1. With this specification, the process has 
jumps of infinite activity (and of bounded variation), which is not covered in the framework of [6]. We 
consider cr = and a = 0.2 so as to study both the bounded and unbounded variation cases. 

We let q = 0.03 and for the surplus cost we consider the quadratic case f{x) = x'^, x E M. By 
straightforward calculation, a = —Cq/2, = —Cq/2 — $(g)~\ 

!^ 2 2s + Cq s^ + Cqs ^. 2 Cq + 2s 

and for x > s 

g*('?)a; 
00 

- J2 B,,,e-^-^^{Cq^'''\s, X- 6,,) + /^^^^ (s, x- e.,,)), 

i=l 

g$((j)x 

^.(^; /') = -'^•(9)) + 2«:(^)(3,a;; -$(g)) 

00 

- B,,,e~^-''-{CqK^^\s, X- e.,,) + 2k'^^\s, x; e^,,)), 
1=1 

where k^") (t, t'; C) := e^^y'^dy for t' > t, C G K and n > 0. 

7.1. Results. For the case K > 0, the first step is to obtain the pair (s*, 5"*) as in Proposition 4.1. As 
is discussed in Section 4.2, starting at s = oq and as we decrease the value of s, we arrive at the desired 
(s*, S*) that makes the function ^(s*, ■) tangent to the x-axis at 5*. Figure 1 plots {G{s, S); S > s} and 
{His, S); S* > s} for s = ao, (ao + s*)/2, s*, -ao/2 + 3s*/2, -oq + 2s*. The lines in red correspond 
to the desired curve; the starting point becomes s* and the point touching zero becomes S*. Here we 
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assume C = K = 10. As it turns out, ^/(s, ■) appears to be strictly convex in this example. Moreover, 
as is already clear analytically, it starts at zero (for the unbounded variation case) or below zero (for 
the bounded variation case). Hence Q{s, S) has a unique global minimum over S E [s, oo), which also 
appears to be increasing in s. Hence, we apply a bisection type method to obtain s* and S*. 




unbounded variation case (cr > 0) 




bounded variation case (a = 0) 

FIGURE 1. Illustrations of {g{s,S);S > s} (left) and {nis,S);S > s} (right) for 

s = ao, (ao + s*)/2, s\ -ao/2 + 3s* /2, -oq + 2s*. 

With (s*, 5"*) computed instantaneously using the technique addressed above, the value function is 
calculated using (5.2). In Figure 2, we first plot it against the initial value x for the proportional cost C = 
30, 20, 10, 5, 1, with the common fixed ordering cost K = 10. The plus signs indicate (s*, Vg^^s* (s*)) 
and {S*, Vs*,s* (S*)) for each C. It can be confirmed that the value function is increasing in C uniformly 
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in X G M. Moreover, s* tends to increase as C decreases. This is consistent with our intuition that one 
is more eager to replenish as the ordering cost decreases. We also see in this plot that S* also tends to 
increase as C decreases. 




unbounded variation case (a > 0) bounded variation case (a = 0) 

Figure 2. The value functions for various values of the proportional cost C. 




unbounded variation case (a > 0) bounded variation case (a = 0) 

Figure 3. The value functions for various values of the fixed cost K. 

We now consider decreasing the value of the fixed cost K and confirm the convergence to the case 
K = as studied in Section 6. In Figure 3, we plot the value functions for = 100, 50, 10, 5, 1 (dotted) 
along with that for the no-fixed cost case given by (6.2) (solid) with the common proportional cost 
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C = 10. The circle signs indicate (ao, ^00(^^0))- We can confirm that, as the value of K decreases, the 
value function converges decreasingly to the one for K = 0. The convergence of the points (s*, S*) to oq 
is also confirmed. Regarding the smoothness of the value function, it appears indeed that it is continuous 
at s* for the case of bounded variation while it is differentiable for the case of unbounded variation. 
Furthermore, the smaller the value of K, the smoother the value function gets. This is consistent with 
Lemma 6.3 where the smoothness holds in a higher order for K = 0. 
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